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ABSTRACT 

We have examined the relationship between rotation and activity in 14 late- type (M6-M7) M dwarfs, 
using high resolution spectra taken at the W.M. Keck Observatory and flux-calibrated spectra from the 
Sloan Digital Sky Survey. Most were selected to be inactive at a spectral type where strong Ha emission 
is quite common. We used the cross-correlation technique to quantify the rotational broadening; six 
of the stars in our sample have usim>3.5 kms -1 . Our most significant and perplexing result is that 
three of these stars do not exhibit Ha emission, despite rotating at velocities where previous work 
has observed strong levels of magnetic field and stellar activity. Our results suggest that rotation 
and activity in late-type M dwarfs may not always be linked, and open several additional possibilities 
including a rotationally-dependent activity threshold, or a possible dependence on stellar parameters 
of the Rossby number at which magnetic/activity "saturation" takes place in fully convective stars. 

Subject headings: stars: low-mass, brown dwarfs — stars: activity — stars: late-type — stars: rotation 



1. INTRODUCTION 

Many M dwarfs, which are the most abundant stars 
in the Milky Way, have strong magnetic dynamos that 
give rise to chromospheric and coronal heating, produc- 
ing emission from the X-ray to the radio. Although 
this magnetic heating (or activity) has been observed 
for decades, the exact mechanisms that control magnetic 
activity in M dwarfs are still not well-understood. In 
the Sun, magnetic fields are generated at the boundary 
between the convective and the radiative zones (known 
as the tachocline), where the differential rotation in the 
convective zone creates a rotational shear (Parker 1993; 
Ossendrijver 2003; Thompson et al. 2003) that allows 
magnetic fields to be generated, stored and ultimately 
rise to the surface. These fields drive the heating of the 
stellar chromosphere and corona, resulting in both flares 
and lower level quiescent magnetic activity. A common 
diagnostic of this heating in optical spectra (of M dwarfs) 
is the Ha emission line. 

Rotation in solar-type stars slows with time due to an- 
gular momentum loss from magnetized stellar winds; as 
a result, magnetic activity decreases. Skumanich (1972) 
found that both activity (as measured by Ca II emis- 
sion) and projected rotational velocity decrease over time 
as a power law (t~ - 5 ). Subsequent studies confirmed 
the Skumanich results and demonstrated that there is a 
strong link between age, rotation and activity in Solar- 
type stars (Barry 1988; Soderblom et al. 1991; Pizzolato 
et al. 2003; Mamajek & Hillenbrand 2008). 

There is strong evidence that the rotation-activity re- 
lation extends from stars more massive than the Sun to 
smaller dwarfs (Pizzolato et al. 2003; Mohanty & Basri 
2003; Kiraga & Stepien 2007). However, at a spectral 
type of ~M3 (0.35 M Q ; Reid & Hawley 2005; Chabrier 
& Baraffe 1997) 4 , stars become fully convective and the 
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tachocline presumably disappears. This transition marks 
an important change in the stellar interior that has been 
thought to affect the production and storage of inter- 
nal magnetic fields. Despite this change in the stellar 
interiors, magnetic activity persists in the late-type M 
dwarfs; the fraction of active stars peaks around a spec- 
tral type of M7 before decreasing into the brown dwarf 
regime (Hawley et al. 1996; Gizis et al. 2000; West et al. 
2004). 

A few previous studies of Ha emission have uncovered 
evidence of a possible rotation-activity relation extend- 
ing past the M3 convective transition and into the brown 
dwarf regime (Delfosse et al. 1998; Mohanty & Basri 
2003; Reiners & Basri 2007). In addition, recent simula- 
tions of magnetic dynamo generation in fully convective 
stars find that rotation may play a role in magnetic field 
generation (Dobler et al. 2006; Browning 2008). How- 
ever, the lack of an unbiased sample of high resolution 
spectra of late-type M dwarfs complicates the situation. 

Using over 30000 spectra from the Sloan Digital Sky 
Survey (SDSS; Adelman-McCarthy et al. 2008), West 
et al. (2006, 2008) showed that the activity fraction of 
M dwarfs varies as a function of stellar age (using Galac- 
tic height as a proxy for age) and that the Ha activity 
lifetime for M6-M7.5 stars is 7-8 Gyr. Nearby samples of 
late-type M dwarfs are therefore biased towards young 
populations with high levels of activity; until recently 
every known M7 dwarf was observed to be magnetically 
active (Hawley et al. 1996; Gizis et al. 2000; West et al. 
2004). 

Due to the intrinsic faintness of late-type M dwarfs, 
previous rotation studies were limited (and therefore 
biased) toward, bright, nearby, and therefore active 
M dwarfs; the 21 M6-M7.5 dwarfs that have previous 
rotation measurements are all active (Delfosse et al. 
1998; Mohanty & Basri 2003; Reiners & Basri 2007). 
The advent of large spectroscopic surveys of M dwarfs 
(e.g. SDSS) allows for improved sample selection that can 

and migrate across a range of spectral types while at the same 
stellar mass; mass and spectral type are only coupled in stars that 
have settled on the main sequence (> 500 Myr for most late- type 
M dwarfs). 
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minimize observational bias in the study of late-type M 
dwarfs. 

In this paper we present results from our study of the 
usini rotation velocities for a small sample of M6-M7 
dwarfs, most of which were selected to be inactive or 
weakly active from the SDSS low-mass star spectroscopic 
sample. We describe the sample and observations in §2, 
our analysis in §3, and the observational results in §4. 
We then discuss our results in §5. 

2. DATA 

Our sample was selected from the West et al. (2008) 
Sloan Digital Sky Survey (SDSS) M dwarf catalog, a 
spectroscopic sample of almost 40000 M and L-type 
dwarfs. We selected the brightest M6 and M7 stars 
which were either inactive (Ha equivalent width < 0.3) 
or weakly active (Ha equivalent width < 1.0). 12 stars 
were selected using these criteria. Two additional active 
(Ha equivalent width > 1.0) M7 dwarfs were added to 
the sample for comparison to previous studies. While our 
sample is not a complete unbiased sample, representative 
of the underlying M dwarf population, it does consist of 
late-type M dwarfs with activity properties selectively 
different than previously observed. 

The 14 stars were observed during four observing runs 
(13 October 2006, 14 October 2006, 25 April 2007, and 
24 January 2008) using the HIRES spectrograph at Keck 
I. Our setup covered the wavelength region from 5600 to 
10000 A with a resolution of R-30000. The data were 
reduced with standard IDL routines that included flat- 
fielding, sky subtraction and cosmic ray removal. In ad- 
dition to our 14 program stars, we observed the nearby 
M6 dwarf G1406, a star with no measurable rotation 
(Delfosse et al. 1998; Mohanty & Basri 2003) that we 
used as a zero rotation velocity template. 

3. ANALYSIS 
3.1. Rotation from Cross Correlation 

To measure the i;sim rotation velocities for our sample, 
we used a cross-correlation technique similar to that of 
previous studies (e.g. Delfosse et al. 1998; Mohanty & 
Basri 2003): we cross-correlated each program spectrum 
with the spectrum of a slowly rotating comparison star. 
The width of the resulting cross correlation function is a 
direct probe of the rotational broadening. 

We used Gl 406 (M6) as our template because it has 
a similar spectral type to our sample as well as a small 
(non detectable) rotation velocity (Delfosse et al. 1998; 
Mohanty & Basri 2003). Four wavelength regions with 
strong molecular bands were used to compute indepen- 
dent cross-correlation functions: 6700-6760 A, 7080-7140 
A, 8430-8500 A (adapted from Browning et al. 2008), 
and 9946-9956 A (Reiners 2007). Figure 1 shows an ex- 
ample of all four regions from one of the stars in our 
sample (SDSS094738.45+371016.5). The S/N ranges 
from 10-60 in all of the spectral regions for all of the 
stars. The three bluer spectral windows yielded good 
cross-correlation functions that were consistent with each 
other for all 14 program stars. In contrast, only a few 
good solutions were produced using the 9946-9956A win- 
dow. The problem with this spectral window may be 
due to the mix of broad and narrow FeH features in the 
9946-9956A region; this region is more amenable to a 



spectral fitting approach to measure rotational broaden- 
ing (e.g. Reiners 2007). Gl 406 is an active M6 dwarf 
with a strong magnetic field (Reiners & Basri 2007), 
and the FeH region is also sensitive to magnetic effects. 
To alleviate this problem, we followed the same cross- 
correlation procedure but instead used one of our non 
rotating (based on the other 3 spectral regions) program 
stars, SDSS125855. 13+052034.7 to cross correlate with 
the sample spectra. Using this template, we were able 
to get better cross-correlation functions for the 9946- 
9956A region. 

To measure the rotational broadening of each spec- 
trum, we compared the resulting cross-correlation func- 
tion to that of a rotationally broadened template. The 
Gl 406 and SDSS125855. 13+052034.7 templates were ro- 
tationally broadened to larger rotation velocities using 
the technique of Gray (1992) and cross-correlated with 
the original (unbroadened) template in 0.5 kms" 1 inter- 
vals. The resulting spun-up cross-correlation functions 
were compared with the cross-correlation functions of 
the sample spectra and a wsini was determined based 
on the best fit spun-up template. Figure 2 shows the 
cross-correlation function of SDSS094738.45+371016.5 
with G1406 in the 7080-7140A region (solid) compared 
with the cross-correlation function of G1406 with the 
best-fit rotationally broadened G1406 spectrum (dotted; 
6 kms" 1 ), and the auto-correlation function of G1406 
(dashed; kms -1 broadening). The cross-correlation 
reveals that SDSS094738.45+371016.5 appears to be ro- 
tating with a velocity > 6 kms -1 . 

No detectable difference was seen between the non- 
rotating template and the spun-up templates with 
V rQ i <3.5 kms -1 for the S/N of our typical spectra. We 
therefore used 4 kms" 1 as our threshold for detectable 
rotation in this paper. The vsim values from each spec- 
tral window were averaged and the resulting value for 
each star is reported in Table 1. For stars with multiple 
spectra, the i>sim values in Table 1 represent the average 
of all spectral windows for all observations. The usini 
variation amongst spectral windows and multiple obser- 
vations was always within the 0.5 kms" 1 bin size. We 
therefore used 0.5 kms" 1 as the uncertainty in the vs'mi 
measurements . 

3.2. Spectral Analysis 
3.2.1. Activity 

All of the spectra were spectral typed by eye using 
the Hammer spectral analysis package (Covey et al. 
2007) on the SDSS spectra. We measured the equiva- 
lent widths (EW) of the Ha emission lines in both the 
SDSS and Keck spectra. The Keck spectra are more 
sensitive to low levels of emission; they can be distin- 
guished better against the pervasive molecular features. 
Our EWs are really psuedo-EWs, since there is no real 
continuum in the spectra. This has different effects at 
low and high resolution, which is one reason that EWs 
tend to be smaller when measured from high resolution 
spectra. Almost all of our targets chosen to be inac- 
tive at low resolution proved inactive even at high res- 
olution, and the EWs when detected were similar. In 
one case (SDSS083231. 52+474807.7), the high resolution 
spectrum seems to have caught a small flare. Figure 3 
shows the SDSS and Keck spectra of the inactive M7 
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Fig. 1 — Keck spectra of SDSS094738.45+371016.5 in the four spectral regions used for cross-correlation: (a) 6700-6760 A, (b) 7080-7140 
A, (c) 8430-8500 A, and (d) 9946-9956 A. 




Fig. 2.— Cross correlation function of SDSS094738.45+371016.5 
with G1406 in the 7080-7140 A region (solid) compared with 
the cross-correlation function of G1406 with the best-fit rotation- 
ally broadened G1406 spectrum (dotted; 6 kms -1 ), and the auto- 
correlation function of G1406 (dashed; kms -1 broadening). The 
cross-correlation reveals that SDSS094738. 45+371016. 5 appears to 
be rotating with a velocity > 6 kms -1 despite not having any signs 
of activity in either the SDSS or Keck spectra. 



dwarf SDSS094738.45+371016.5 (a and b respectively) 
and the active M7 dwarf SDSS162718. 20+353835. 7 (c 
and d). Dashed lines in the SDSS spectra (a and c) 
indicate the expected position of the Ha emission line. 
The inset panels in the SDSS spectra are zoomed in on 



the Ha region (6500-6600 A). Despite their high quality, 
it is clear that no Ha emission is detected in either of the 
inactive spectra. 

We also measured the H/3, H7, R5, and Call K emis- 
sion lines from the SDSS spectra. The EW values are 
presented in Table 1. For the stars in which no lines 
were detected, we computed EW upper limits based on 
the required EW needed to detect an emission line at 
the 3cr confidence level. These are also reported in Table 
1. Using the relations from Walkowicz et al. (2004, Ha) 
and West & Hawlcy (2008, other lines), we converted 
the EWs and upper limits for all of the lines to fractions 
of the bolometric luminosity (e.g. Lna/Lboi, etc.). The 
Lii ne /Lboi values are presented in Table 2. 

We estimated the Rossby number (i? =Period/r c ) for 
the 6 stars with measured rotation velocities, following 
the prescription of Reiners et al. (2008, r c =70 days), and 
assuming stellar radii of 0.12 R Q and 0.1 R© for M6 and 
M7 dwarfs respectivley (Reid & Hawley 2005). Our es- 
timated Rq values are included in Table 1. All of the 
estimated Rossby numbers are in the regime of activ- 
ity "saturation" (R < 0.1) described in Reiners et al. 
(2008). 

3.2.2. Metallicity 

We measured the Ti05, CaH2, and CaH3 molecular 
band indices (Reid et al. 1995) from the SDSS spec- 
tra, features that probe the relative metallicities of M 
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TABLE 1 

Measured Attributes 



TVT 

IN 3,1X16 










Equivalent Width (A) 








1 1 wo 


v^axiz 




rtossby 






(kms -1 ) 


Hct a 




Ha b 


H/3 a 


±± j 


H<5 a 




Call K a 










Num. 


SDSS011012. 22-085627. 5 


M7 


<3.5 


<0 


2 


<0.2 


<0.6 


<1.7 


<4 


3 


<0 


9 


0.17 


0.27 


0.56 




SDSS021749. 99-084409. 4 


M6 


4.0+0.5 


<0 


1 


<0.1 


<0.4 


<2.2 


<1 


7 


<0 


8 


0.28 


0.35 


0.67 


0.03 


SDSS023908. 41-072429. 3 


M7 


<3.5 


<0 


1 


0.3 


<0.8 


<7.0 










0.23 


0.32 


0.61 




SDSS072543.94+382511.4 


M7 


<3.5 


<0 


1 


<0.1 


<0.2 


<0.5 


<0 


5 


<1 


1 


0.27 


0.35 


0.67 




SDSS083231. 52+474807.7 


M6 


4.0±0.5 


<0 


1 


3.9 d 


<0.2 


<0.5 


<0 


5 


1 


1 


0.28 


0.34 


0.64 


0.03 


SDSS094720.07-002009.5 


M7 


6.5±0.5 


3 


4 


4.0 


4.2 


17.2 


2 


6 


12 


2 


0.22 


0.28 


0.57 


0.01 


SDSS094738.45+371016.5 


M7 


6.0+0.5 


<0 


2 


<0.1 


<0.2 


<0.9 


<1 


8 


<1 





0.19 


0.27 


0.59 


0.01 


SDSS110153.86+341017.1 


M7 


<3.5 


<0 


2 


<0.1 


<0.2 


<1.1 


<1 


1 


<0 


8 


0.22 


0.29 


0.60 




SDSS112036.08+072012.7 


M7 


<3.5 


<0 


1 


<0.1 


<0.2 


<0.3 


<0 


4 


<0 


5 


0.26 


0.32 


0.66 




SDSS125855.13+052034.7 


M7 


<3.5 





9 


<0.1 


1.8 


7.2 


4 


5 


5 


8 


0.25 


0.31 


0.62 




SDSS151727.72+335702.4 


M7 


4.5±0.5 


<0 


2 


<0.1 


<0.2 


<0.6 


<0 


8 


<0 


3 


0.26 


0.35 


0.70 


0.02 


SDSS162718.20+353835.7 


M7 


8.0+0.5 


5 


9 


5.4 


19.7 


37.1 


32 


4 


26 


5 


0.24 


0.20 


0.61 


0.01 


SDSS220334.10+130839.8 


M7 


<3.5 


<0 


1 


<0.1 


<0.6 


<2.8 


<2 


3 


<2 


3 


0.22 


0.32 


0.66 




SDSS225228.50-101910.9 


M7 


<3.5 


<0 


2 


<0.2 


<0.7 


<9.8 










0.19 


0.26 


0.59 





a Equivalcnt width or 3c upper limit measured from the SDSS spectrum 
^Equivalent width or 3(7 upper limit measured from the Keck spectrum 
c Bandhcad measured from the SDSS spectrum 
d Keck spectrum observerd during a small flare 
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Fig. 3.— SDSS and Keck spectra of the inactive M7 dwarf SDSS094738.45+371016.5 (a and b respectively) and the active M7 dwarf 
SDSS162718. 20+353835. 7 (c and d). Dashed lines in the SDSS spectra (a and c) indicate the expected position of the H« emission line. 
The inset panels in the SDSS spectra are zoomed in on the Ho region (6500-6600 A). Despite their high quality, it is clear that no Ha 
emission is detected in either of the inactive spectra. 



dwarfs (Lepine et al. 2003; Burgasser & Kirkpatrick 2006; 
West et al. 2008). The combination of these indices 
(CaH2+CaH3 vs. Ti05) has been used to separate the 
dwarf population ([Fe/H]~0) from the subdwarf popula- 
tion ([Fe/H]~-1). We include the index value for each 
molecular band in Table 1. All of the index values are 
consistent with disk dwarf metallicity (near Solar). 



3.2.3. Dynamics 

We assigned a Galactic population to each star using 
the 3-D space motions computed in West et al. (2008) 
and the criteria of Leggett (1992), where young disk stars 
(YD) have -20 < U < 50, -30 < V < 0, and -25 < W < 
10. Stars outside this ellipse or with |W| > 50 arc old 
disk (OD) and stars near the edge but with |W| < 50 are 
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classified as young-old disk stars (YOD). All of the dy- 
namic classifications are included in Table 1. While our 
sample includes all 3 of the components, it is dominated 
by the old disk and yound old disk populations; 12 of 14 
stars are OD or YOD. 

4. OBSERVATIONAL RESULTS 

Six of the fourteen M6-M7 dwarfs in our sample have 
detectable rotation. Three of the rotating stars have 
measurable activity but the other 3 show no signs of ac- 
tivity in any of the emission lines in cither the SDSS or 
Keck spectra. The cross-correlation shown in Figure 2 
(SDSS094738.45+371016.5) is an example of one of the 
inactive M7 dwarfs that appears to be rotating despite 
not being magnetically active. 

Figure 4 shows Lua/Lboi (activity) as a function of 
vsini for the M6-M7.5 dwarfs from this paper (filled sym- 
bols) and previous studies (Delfosse et al. 1998; Mohanty 
& Basri 2003; Reiners 2007, open diamonds). Lna/Lboi 
values were calculated from EW measurements using the 
X conversions of Walkowicz et al. (2004). Our sam- 
ple includes M6-M7 dwarfs with both measured rota- 
tion as well as activity from the SDSS and Keck spec- 
tra (filled circles), measured rotation and activity from 
the Keck spectra (Call K detected in SDSS; filled dia- 
monds), activity from the Keck spectra but no rotation 
(filled triangles), activity from the SDSS spectra but no 
rotation (filled squares) , no rotation or activity (horizon- 
tally aligned dots) and measured rotation but no activity 
(filled stars). Upper limits in both velocity and activ- 
ity denote the levels to which our sample (and previous 
studies) could probe. All previous M6-M7.5 dwarfs were 
found to be active, while 9 of the 14 stars in our sample 
show no activity in either the SDSS or Keck spectra. The 
lack of activity is not surprising since that was our main 
selection criterion, however it highlights the fact that we 
are probing a sample with very different properties than 
previously studied. 

The two most active stars in our sample are also the 
fastest rotators, and 6 of the inactive stars have no de- 
tectable rotation. Two of the active stars show no sign 
of rotation, which may be due to inclination effects (al- 
though the incidence in that case would be higher than 
statistically expected, if such small numbers can be said 
to have any significance). 5 A similar excess of active non- 
rotators was found by Browning et al. (2008) at spectral 
types M3-M5. All 6 of the stars with detectable rotation 
have estimated Rossby numbers that are beyond (smaller 
than, by an order of magnitude) the activity "saturation" 
threshold of i?o=0.1 discussed in Reiners et al. (2008, 
the stars would need to have Y ro t <0.7 kms -1 to have 
Ro >0.1). However, none of these rotators show activ- 
ity at the "saturated" level of Lna/Lboi ~ 10~ 4 . Most 
intriguing are the three M7 dwarfs that show no activ- 
ity but are clearly detected rotators. We discuss this in 
more detail in §5. 

Neither the metallicity nor dynamic population analy- 
ses reveal any significant trends with either rotation ve- 
locity or activity. This is not surprising, since both mea- 

5 Using the statistical formalism outlined in Browning et al. 
(2008), the probability of selecting 2 out of 5 active stars with the 
inclination required to alter their observed vsini velocities from the 
mean detected velocity of 6 kms -1 to 3.5 kms -1 (non-detection) 
is 24%. 



surements are only relevant in a statistical sense. How- 
ever, all but one of the stars with a detected rotation 
comes from the YD or YOD population, and the one ro- 
tating OD shows no signs of activity (as do one of each 
from the other two populations). 

5. DISCUSSION 

We conducted high resolution spectral observations of 
14 M6-M7 dwarfs and found vsini rotation velocities for 
6 for the stars. Three of the stars showed both activity 
and rotation, 6 of the stars showed neither rotation nor 
activity, 2 of the stars showed activity but no rotation 
and 3 stars showed rotation but no activity. These results 
are in contrast with previous studies that found a strong 
connection between rotation and activity in all (active) 
M6-M7 dwarfs (Delfosse et al. 1998; Mohanty & Basri 
2003; Reiners 2007). As can be seen in Figure 4, those 
studies produced one out of 20 stars which was active but 
not rotationally broadened (very reasonably explained by 
the inclination effect; there is a 94% chance of selecting 
one out of 20 with the inclination needed to take it from 
4 kms -1 to 2 kms -1 ), and one which showed rotation 
and low (but not zero) activity. Our sample is the first 
rotation study to include by design M6-M7 dwarfs that 
are inactive. 

Whereas previous studies sampled a range of dynamic 
populations and therefore stellar ages (Delfosse et al. 
1998; Mohanty & Basri 2003), our sample was selected 
from a more distant, mostly inactive and potentially 
older population (M7 dwarfs have activity lifetimes of 
~8 Gyr; West et al. 2008). Although our sample is not 
representative of the entire M6-M7 population, it does 
examine a different part of parameter space than pre- 
viously sampled and observes something not previously 
seen: rotation in inactive M7 dwarfs. Currently, there 
is no clear explanation for the rotation-activity patterns 
seen in our sample. However, we present a number of 
possible explanations that may account for this unusual 
dataset. 

One possible explanation is that rotation and activity 
are not strongly linked in the fully convective envelopes 
of M6-M7 dwarfs. Previous studies only examined active 
M6-M7 dwarfs, many of which have activity levels with 
no notable correlation with rotation beyond a general 
tendancy for rotators to lie in a broadly spread "sat- 
urated" activity regime. In fact, several active dwarfs 
have no detectable rotation (see Figure 4). Besides satu- 
ration, another explanation for the lack of a strong rela- 
tionship between activity and rotation is inclination; the 
rotation axes of the stars may be inclined and therefore 
have reduced line-of-sight components to their rotation 
velocities. On the other hand, the 3 inactive stars with 
detected rotations certainly cannot be resolved by sim- 
ply imposing inclination effects, which can only serve to 
make their true rotations even higher. 

Another possibility is that there is a rotationally- 
induced threshold for activity at V ro t >6 kms -1 (Haw- 
ley et al. 1996). This would explain the inactive rotators 
(they are below the threshold) and require that active 
stars with vsini rotation velocities less than the thresh- 
old have inclined rotation axes. Combining all of the 
previous rotation studies for M6-M7 (Figure 4), we see 
that 7 of the 26 active dwarfs have vsini < 6 kms -1 ; the 
rest of the stars are rotating several km s _1 above the de- 
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TABLE 2 
Measured Attributes II 



Name 


Spec. 
Type 


vsini 
(kms- 1 ) 


LHa/L6o! a 

(xlO- 6 ) 


LHa/Lf,o! b 

(xlO- 6 ) 


LH/3/Lto! a 

(xlO- 6 ) 


(xlO- 6 ) 


LH«/Lf,o! a 

(xlO- 6 ) 


LcaIIK/L6o! a 

(xlO- 6 ) 


Dyn. c 
Pop. 


SDSS011012.22-085627.5 


M7 


<3.5 


<1.0 


<1.0 


<0.5 


<0.3 


<0.9 


<0.2 


YOD 


SDSS021749.99-084409.4 


M6 


4.0±0.5 


<2.0 


<2.0 


<1.4 


<2.4 


<1.5 


<0.8 


YD 


SDSS023908.41-072429.3 


M7 


<3.5 


<0.5 


1.6 


<0.6 


<1.4 






YOD 


SDSS072543.94+382511.4 


M7 


<3.5 


<0.5 


<0.5 


<0.2 


<0.1 


<0.1 


<0.2 


OD 


SDSS083231. 52+474807.7 


M6 


4.0±0.5 


<2.0 


69.0 


<0.7 


<0.6 


<0.5 


1.1 


OD 


SDSS094720.07-002009.5 


M7 


6.5±0.5 


18.0 


21.0 


3.4 


3.4 


0.52 


2.4 


YOD 


SDSS094738.45+371016.5 


M7 


6.0±0.5 


<1.0 


<0.5 


<0.2 


<0.2 


<0.4 


<0.2 


OD 


SDSS110153.86+341017.1 


M7 


<3.5 


<1.0 


<0.5 


<0.2 


<0.2 


<0.2 


<0.2 


OD 


SDSS112036.08+072012.7 


M7 


<3.5 


<0.5 


<0.5 


<0.2 


<0.1 


<0.1 


<0.1 


OD 


SDSS125855.13+052034.7 


M7 


<3.5 


4.68 


<0.5 


1.4 


1.4 


0.90 


1.2 


OD 


SDSS151727.72+335702.4 


M7 


4.5±0.5 


<1.0 


<0.5 


<0.2 


<0.1 


<0.2 


<0.1 


YOD 


SDSS162718.20+353835.7 


M7 


8.0±0.5 


31.0 


28.0 


16.0 


7.4 


6.5 


5.3 


YOD 


SDSS220334.10+130839.8 


M7 


<3.5 


<0.5 


<0.5 


<0.5 


<0.6 


<0.5 


<0.4 


YD 


SDSS225228.50-101910.9 


M7 


<3.5 


<1.0 


<1.0 


<0.6 


<2.0 






YOD 



NOTE. — The L^ne/Lbd values were calculated using the \ conversion factors derived by Walkowicz ct al. (2004) and West & Hawley (2008). 

a Equivalcnt width or 3<r upper limit measured from the SDSS spectrum 
b Equivalent width or 3cr upper limit measured from the Keck spectrum 
c Dynamical Populations based on the criteria in Lcggctt (1992) 
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Fig. 4. — LHa/Lbol (activity) as a function of vs'mi (rotation) for the M6-M7.5 dwarfs from this paper (filled symbols) and previous 
studies (Dclfossc et al. 1998; Mohanty & Basri 2003; Rcincrs 2007, open diamonds). hn a /h bo i values were calculated from equivalent width 
measurements using the \ conversions of Walkowicz et al. (2004). Our sample includes M6-M7 dwarfs with both measured rotation as 
well as activity from the SDSS and Keck spectra (filled circles), measured rotation and activity from the Keck spectra (Call K detected in 
SDSS; filled diamonds), activity from the Keck spectra but no rotation (filled triangles), activity from the SDSS spectra but no rotation 
(filled squares), no rotation or activity (horizontally aligned dots) and measured rotation but no activity (filled stars). Upper limits in both 
velocity and activity denote the levels to which our sample (and previous studies) could probe. All previous M6-M7.5 dwarfs were found 
to be active, while 9 of the 14 stars in our sample show no activity in either the SDSS or Keck spectra. The dearth of inactive late-type M 
dwarfs in previous studies is due to a selection effect that biases nearby samples to younger, more active stars (West et al. 2006, 2008). 3 
of our stars show strong evidence for rotation despite having no activity. Their usim velocities are well above the 0.7 kms -1 required to 
have Ro > 0.1 and be beyond the "saturated" regime (Reiners et al. 2008). 
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tection thresholds. To quantify this possibility we follow 
the statistical analysis of Browning et al. (2008), where 
we compute the probability (from the binomial distribu- 
tion) of observing n objects with inclination i < i cr it- We 
find that the probability of randomly drawing a sample 
of 7 stars with the inclinations necessary to take the slow 
but active stars from a V ro t=6.5 kms -1 to the observed 
vsvai is 85.7%. If we instead calculate the more conser- 
vative probability of drawing inclinations that take the 7 
stars from 8 kms' 1 (the median vsim) to their respective 
vsvai values, we obtain 55.9%. In this scenario, the inac- 
tive stars, which are presumably older, have spun down 
enough to fall below the rotation-activity threshold. 

In earlier M dwarfs, the mere detection of rotation im- 
plies that the star is in the "saturation" regime; an ef- 
fect that should have become even more pronounced for 
these even smaller objects (assuming the convective over- 
turn time is similar). All of our detected rotators have 
estimated Rossby numbers that are beyond the "satura- 
tion" threshold for early-type M dwarfs (Reiners et al. 
2008). One possibility is that in these fully convec- 
tive late-type M dwarfs, the Rossby number "saturation" 
threshold might be substantially smaller. However, this 
model would imply that stars go from fully active to com- 
pletely inactive over a small range of Rossby numbers - 
something that contradicts what is seen in early-type M 
dwarfs, where the activity falls smoothly with increasing 
Rossby number (Kiraga & Stepien 2007; Reiners et al. 
2008). It is clear that the velocity thresholds proposed 
here would have to be lowered (as a function of spectral 
type) to be compatible with the results of Reiners et al. 
(2008); too many of their stars are active below vsim = 
6 kms -1 (including several in the M5-6 spectral range). 

Of course one additional possibility is that M6-M7 
dwarfs might have activity cycles and that we are ob- 
serving the inactive stars at a minimum in their cycles. 
Although this might be observed for a small number of 
stars, it is unlikely that an activity minimum can explain 
the majority of inactive M dwarfs seen in the Milky Way. 
West et al. (2008) showed that the activity fractions of M 
dwarfs are strongly correlated with the height above the 
Galactic plane (and therefore age), a trend that cannont 
be explained by an activity cycle. 

The dearth of a large sample of nearby, inactive late- 
type dwarfs precludes large-scale follow-up that would 
further test the age-rotation-activity relation of fully con- 
vective low-mass dwarfs. However, future efforts to ex- 
tend our analysis to both inactive M5 and M8 dwarfs 
will help elucidate the problem. The multi-epoch pho- 
tometric observations produced by all sky time-resolved 
surveys (e.g. Pan-STARRS and LSST) will produce ro- 
tation periods for a large number of late-type M dwarfs 
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that will eliminate the inclination effects. Of course, in- 
active M dwarfs present a problem to this method in 
that they may not exhibit photometric variability (if they 
don't produce starspots) . The combination of these and 
other future observations will help advance the under- 
standing of what role rotation plays in the production of 
magnetic fields in fully convective objects, and provide 
crucial constraints for theoretical models of fully convec- 
tive dynamos. They will also be a key component to 
understanding the feedback between magnetic fields and 
stellar angular momentum evolution through stellar ac- 
tivity. 
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